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We extend simulations of simplicial gravity in four dimensions to include degenerate triangulations and demon- 
strate that using this ensemble the geometric finite-size effects are much reduced. We provide strong numerical 
evidence for the existence of an exponential bound on the entropy of the model and establish that the phase 
structure is identical to that of a corresponding model restricted to an ensemble of combinatorial triangulations. 



1. INTRODUCTION 

In discretized models of four-dimensional Eu- 
clidean quantum gravity, known as simplicial 
gravity, the integration over metrics is replaced 
by summations over an ensemble of triangulations 
constructed by all possible gluings of equilateral 
4-simplexes into closed (piece-wise linear) simpli- 
cial manifolds (see e.g. Ref. Q). The regularized 
Euclidean Einstein-Hilbert action is particularly 
simple in this approach; it can be taken to de- 
pend on only two coupling constants, k and n, 
related to the inverse Newton's and cosmological 
constants. The regularized grand-canonical par- 
tition function thus becomes: 
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The sum is over all distinct triangulations TeT, 
Nt is the number of i-simplexes in the triangula- 
tion T and Ct denotes its symmetry factor — the 
number of equivalent labelings of the vertexes. 

Extensive numerical simulations have estab- 
lished that the model Eq. (|l|) has a strong- 
coupling (small k) crumpled phase and a weak- 
coupling (large k) elongated phase, separated by 
a discontinuous phase transition. In the crum- 
pled phase the geometry is dominated by a sin- 
gular structure; two singular vertexes connected 
to an extensive fraction of the total volume, 
joined by a sub-singular edge. The elongated 
phase, on the other hand, is dominated by essen- 
tially one-dimensional (tree-like) triangulations 
— branched polymers. 



2. DEGENERATE TRIANGULATIONS 

In Eq. ([!]), T denotes a suitable ensemble of 
triangulations included in the partition function. 
Different ensembles arc defined by imposing var- 
ious restriction on how the simplexes are glued 
together. Provided this leads to a well-defined 
partition function, and as long as this difference 
is only at the level of discretization, one expects 
different choices of T to lead to the same contin- 
uum theory in the thermodynamic limit. This is 
known to be true in two dimensions where mod- 
els of simplicial gravity corresponding to different 
choices of T are soluble as matrix models j| . 

All simulations of four-dimensional simplicial 
gravity have, as of today, used an ensemble of 
combinatorial triangulations Tq . In a combinato- 
rial triangulation every 4-simplex is uniquely de- 
fined by a set of five distinct vertexes — it is said 
to be combinatorially unique. Here we report on 
simulations of the model Eq. ([j]) defined with a 
larger ensemble including degenerate triangu- 
lations. We relax the above constraint and allow 
distinct simplexes to be defined by the same set 
of vertexes. We do, however, retain the restric- 
tion that every 4-simplex is defined by a set of 
distinct vertexes, i.e. we exclude degenerate sim- 
plexes. Clearly Tq C T-p. 

The benefit of using a larger ensemble of trian- 
gulations is well known from simulations of sim- 
plicial gravity in two |3| and three dimensions. 
It is established that less constrained the triangu- 
lations are the smaller the geometric finite-size ef- 
fects are. Simulations of four-dimensional simpli- 
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cial gravity are notoriously time-consuming, pri- 
marily due to the large volumes needed to ob- 
serve any "true" infinite volume behavior, hence 
any reduction in the finite-size effects is of great 
practical importance. 

3. RESULTS 

We have simulated the model Eq. (|l|) using de- 
generate triangulations on volumes up to 25.600 
4-simplexes using Monte Carlo methods. As cus- 
tomary we work in a quasi-canonical ensemble of 
spherical manifolds with almost fixed N4: 



Z(p, k; JV 4 ) 
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where Qn 4 (k) = J2TeT(N 4 ) ex P( K ^) ^ s the 
canonical partition function. As there do not ex- 
ist ergodic volume conserving geometric moves, 
the canonical ensemble cannot be simulated di- 
rectly, we must allow the volume to fluctuate. 
The quadratic potential term added to the ac- 
tion ensures, for an appropriate choice of <5, that 
these fluctuations are small. 

In the simulations the triangulation-space is 
explored using a set of local geometric changes, 
the (p, g)-moves [gj. For combinatorial trian- 
gulations the (p, g)-moves are known to be er- 
godic for D < 4. To demonstrate that the same 
holds true for degenerate triangulations we ob- 
serve that every set of combinatorially equivalent 
simplexes, or sub-simplexes, can be made distinct 
by a finite-sequence of the (p, q)-moves. Thus ev- 
ery degenerate triangulations can be reached from 
a combinatorial one. In addition, the local na- 
ture of the (p, g)-moves prohibits the creation of 
pseudo-manifolds in the simulations, i.e. triangu- 
lations containing vertexes with a neighborhood 
not homeomorphic to the 4-sphere. 

A major benefit of including degenerate trian- 
gulations is the reduction of geometric finite-size 
effects. An example of this is the volume depen- 
dence of the pseudo-critical cosmological constant 
Mc(A4) shown in Figure 1. For comparison we 
show the corresponding values for combinatorial 
triangulations. For degenerate triangulations we 
observe a rapid convergence to an infinite volume 
value p; this is quantified by fits to two different 
functional forms: a power-law convergence, 




Figure 1. The pseudo-critical cosmological con- 
stant, p^(N^), together with fits to power-law 
convergence (solid) and to logarithmic divergence 
(dashed). Included are the corresponding values, 
Pc{Na), for combinatorial triangulations. 

Table I 

Parameters in fits of p®{N±) to Eqs. (||) and (0), 
including volumes N4 = 400 to 25.600. 



p 



7 



Eq. 
Eq. 



X7(d.0.f). 



(3) 2.556(3) 0.55(5) 3.8 

(4) 2.385(4) 117 



/i c (iV 4 ) = p+ J^f, 

and a logarithmic divergence, 
p c (N 4 ) = p. + b\ogN 4 . 



(3) 



(4) 

The fit parameters and the quality of the fits are 
shown in Table 1. Contrary to combinatorial tri- 
angulations, for degenerate triangulations there is 
no comparison in the quality of the fits; the latter, 
corresponding to a super-exponential growth of 
the canonical partition function, is ruled out by a 
X 2 /(d.o.f.) w 117 for a fit to Eq. (|). In contrast, 
for combinatorial triangulations it is difficult to 
distinguish between the two scenarios based on 
the quality of the fits [||. The importance of 
this result is that it provides strong numerical 
evidence for an exponential bound on the ensem- 
ble of degenerate triangulations as a function of 
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Figure 2. The (normalized) distributions of ver- 
texe orders for degenerate triangulations. This is 
for N4 > 400 and k = 0. Insert: The number of 
singular vertexes n s . 



volume. And, as Tq 6 7p, this implies an ex- 
ponential bound on the number of combinatorial 
triangulations as well. 

As for combinatorial triangulations, in the 
crumpled phase the internal geometry of degen- 
erate triangulations is dominated by a singular 
structure. The probability distribution of the ver- 
tex orders Oj — the number of 4-simplexes con- 
taining the vertex — contains an isolated peak in 
the tail, indicating singular vertexes (Figure 2). 
However, the distributions P{o) differs in two re- 
spects from the corresponding ones measured on 
combinatorial triangulations: 

(i) The number of singular vertexes is larger 
than two and increases (logarithmically) with 
the volume, i.e. a gas of singular vertexes. 

(ii) On each volume the distribution effectively 
separates in two, depending on whether the 
vertex orders are odd or even. 

It is not clear though, how much significance 
should be attached to this observed difference. 
Due to the collapsed nature of the internal ge- 
ometry it is unlikely that any sensible continuum 
limit exist in the crumpled phase, hence there is 
no reason to expect identical scaling behavior for 
the two different ensembles. 

Additional evidence of a collapsed intrinsic ge- 
ometry comes from the (absence of) volume scal- 



ing of the simplex-simplex correlation function, 
from which we conclude that du = 00. 

We have also investigated the phase structure 
of the model for non-zero values of the inverse 
Newtons's constant k. As for combinatorial trian- 
gulations we observe a phase transition at a value 
k c w 1.5. For k > n c the model is in a branched 
polymer phase; this we establish by measuring the 
fractal dimensions dn and the spectral dimension 
d s for k — 2.0. Including measurements on vol- 
umes, N 4 = 400 to 1600 we get d H = 1.9(2) and 
d s = 1.32(5), in excellent agreement with d# = 2 
and d s =4/3 expected for branched polymers. 

Our results demonstrate that including de- 
generate triangulations in simulations of four- 
dimensional simplicial gravity has many potential 
advantages over a model restricted to combinato- 
rial triangulations. This agrees with the same ob- 
servation previously made in both two and three 
dimensions. The chief benefit is the reduction 
in geometric finite-size effects, mainly due to an 
enlarged ensemble — with a larger triangulation- 
space the infinite- volume fractal structure is more 
easily approximated on small volumes. Further 
work is needed to fully explore the nature of this 
ensemble, for now the most important result we 
want to emphasize is the strong numerical evi- 
dence for an exponential bound on the canonical 
ensemble f2jv 4 of degenerate triangulations. 
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